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I. Introduction 



In this paper wc consider the Anisotropic Mancv Problem (AMP) that was 
introduced by Diacu^ in the early 1990s. The work on the AMP was inspired 
by the Anisotropic Kepler Problem introduced by Gutzwiller in the early 
1970s. Gutzwiller aimed to find connections between classical and quan- 
tum mechanics. His interest was stimulated by an old unsolved quantum 
mechanical problem formulated in a paper written by Einstein:^ even if the 
Born-Sommerfeld-Einstein condition (e.g. see Ref. 2) were appropriate to 
describe the semi-classical limit of quantum theory it was unclear how to find 
a classical approximation for nonintegrable systems. 

Similarly the main reason for considering the AMP is to further analyze 
similarities between classical mechanics and quantum theory. Moreover, as it 
was remarked in Ref. 1, the AMP also brings general relativity into the game, 
since, the Manev's potential, explains the perihelion advance of the inner 
planets with the same accuracy as general relativity.^ It should be remarked 
that bringing general relativity into the game is of particular importance 
since a satisfactory quantum theory of gravitation does not exist. 
Some of the qualitative features of the Anisotropic Manev Problem have 
already been studied. In Ref. 1, a large class of capture-collision and ejection- 
escape solutions is studied by means of the collision and infinity manifold 
techniques. In particular that paper also brought arguments favouring the 
chaoticity and nonintegrability of the system by showing the existence of 
heteroclinic orbits within the zero energy manifold. In Ref. 4 the occurrence 
of chaos on the zero energy manifold and the nonintegrability are finally 
proved, putting into evidence that the AMP is a very complex problem. 
In this work, to gain a better understanding of the complicated dynamics of 
the AMP, we find the symmetric periodic orbits. Analyzing those orbits is 
especially important since, by now, it is well known that studying periodic 
orbits is a valuable general approach to tackle complex problems in classical 
mechanics. The existence of periodic orbits for small values of the anisotropy 
is proved using generalizations of the Poincare continuation method devel- 
oped in Refs. 5-7. 
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The (planar) anisotropic Manev problem is described by the Hamiltonian 
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(1) 



where /i > 1 is a constant and q = (x, y) is the position of one body with 
respect to the other considered fixed at the origin of the coordinate system, 
and p = {Px,Py) is the momentum of the moving particle. The constant 
measures the strength of the anisotropy and for /x = 1 we recover the classical 
Manev problem. Furthermore the equation of motion can be expressed as 



Now consider weak anisotropies, i.e choose the parameter jj, > 1 close to 1. 
Introducing the notation r = a/x^ + and e = /i — 1 with e -C 1 we can 
expand the Hamiltonian (|^) in powers of e and obtain 



It should be pointed out that the term W{r, 9) becomes unbounded as r — >■ 
so that a perturbation analysis is not correct on the ejection-collision orbits. 
This means that the global dynamics of the AMP cannot be completely 
described by perturbations to the Manev problem even at the limit e — 0. 
However many interesting results concerning the Hamiltonian ([^) for weak 
anisotropies (i.e. e -C 1) can be found studying the Hamiltonian (^, some 
of which are presented in this paper. 

In the next section we describe the symmetries of the AMP and we find some 
properties that will be useful to find symmetric periodic orbits. In Section HI 
we prove a continuation theorem for the symmetric periodic orbits of "second 
kind", i.e. the non-circular ones. In Section IV we prove a continuation 
theorem for the orbits of "first kind", i.e. the circular ones, following the 
method developed in Ref. 8. 

II. Symmetries of the Anisotropic Manev Problem 

To find periodic orbits in the anisotropic problem it is peculiarly important to 
know the symmetries of the system, as it was, for example observed in Refs. 




(2) 



H 




cos^9 = Ho + eW{r,e). (3) 
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5, 6. The symmetries of the problem under discussion have been examined in 
Ref. 1 and they are the same as the ones found in Ref. 7 for the anisotropic 
Kepler problem: 



E 

So 
Si 

S2 

S3 
S4 
S5 
Sh 



{x,y,Px,Py,t) 

ix,y,Px,Py,t) 

(x, y, Pxi Py, t) 

(x, y, Pxi Py, i) 
(x, ?/, Pxi Pyi i) 

{XiViPxiPyit) 
(XiViPxiPyit) 



(XiyiPxiPyit) 

{xiyi -pxi -Pyi -t) 
(x, -y, -PxiPyi -t) 
{-XiyiPxi-Pyi-t) 

{-x,-yi -Pxi-Pyit) 

{-X,y, -PxiPyit) 
(X, -ViPxi -Pyit) 

(-X, -yi Pxi Pyi -t) 



(4) 



where E is the identity. 

The symmetries above can be interpreted in the following way: let 7(t) be 
a solution of (H), then S'j(7(t))is another solution for i G {0,1,2,3,4,5,6}. 
For i G {0, 1, 2, 3, 4, 5, 6} the orbit 7(t) will be called symmetric if and only 
if5,(7(t))=7W- 

Let us remark that the symmetries in (^), together with the composition of 
functions, denoted by o, form an abelian group in which the operation acts 
according to the table below. 
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From the table above it is easy to deduce the following 

Proposition 1 The symmetries of the Anisotropic Manev Problem form an 
elementary abelian group of order eight, i.e. a group isomorphic to Z2 x Z2 x 
Z2. 

The symmetries in (H), (except E and 5*6) are very useful to find symmetric 
periodic orbits, especially by means of the continuation method, as we show 
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in the next two sections. Some important properties of the symmetric orbit, 
summarized in Ref. 7 , are expressed in the following lemma: 

Lemma 1 (i) For i — 1 (resp. i = 2) we have that an orbit 7(t) is Si- 
symmetric if and only if it crosses the x axis (resp. y axis) orthogonally. 

(a) An orbit ^{t) is So-symmetric if and only if it has a point on the zero 
velocity curve. 

(Hi) For i = 4,5 an orbit ^{t) is Si-symmetric if and only if it is So- 
symmetric. 

(iv) All the Sz-symmetric periodic orbits are periodic. 

The properties of the y^j-symmetric orbits were first studied by Birkhoff^ 
for the restricted three body problem and later by many other authors. In 
particular Casasayas and Llibre'' state a proposition that gives a technique 
useful to obtain symmetric periodic orbits with respect to 5*0, 5*1, 5*2 for the 
anisotropic Kepler problem that are verified also for the problem under dis- 
cussion in this paper: 

Proposition 2 (i) Fori = 1 (resp. i = 2) we have that an orbit ^(t) is a 
Si-symmetric periodic orbit if and only if it crosses the x axis (resp. y 
axis) orthogonally at two distinct points. 

(a) An orbit 7(t) is a So-symmetric periodic orbit if and only if it meets 
the zero velocity curves at two distinct points. 

(Hi) An orbit j(t) is a Si and S2 symmetric periodic orbit if and only if it 
crosses the x axis and the y axis orthogonally. 

(iv) Fori = 1,2 an orbit ^{t) is a So and Si-symmetric periodic orbit if and 
only if it meets the zero velocity curve and crosses the x, respectively y 
axis orthogonally. 

(v) For i = 4,5, if an orbit 7(t) is Si-symmetric then it is So-symmetric 
and periodic. 

Now we want to find the symmetric periodic orbit for the unperturbed 
problem (e = or = 1) and continue them to periodic solutions of the 
anisotropic system (for e <^ 1). Firstly we observe that, by Proposition 2, 
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the Si symmetric orbits with i — 0, 4, 5 must meet the zero velocity curve 
at two points, i.e. there must be a point where K = 0, but since i^' is a 
constant of motion it must be zero along the orbit. Therefore such orbits 
are ejection-collision orbits, are not periodic and cannot be studied by means 
of the continuation method. Hence we are going to consider the symmetric 
periodic orbit with i = 1,2, and also the ones with i — 3 that are the circular 
orbits of the unperturbed problem. 

To exploit those properties of the symmetric periodic orbits it is convenient 
to write the equation of motion in different coordinates. 
For the Si symmetric orbits with i = 1, 2, as it was noted in Ref. 5 , it is 
convenient to write the canonical equations of the restricted three body prob- 
lem using the Delaunay variables in the rotating frame. ^ Also the Poincare 
synodic variables can be used to find symmetric periodic orbits of the re- 
stricted three body problem.^ The anisotropic Manev problem is different 
since the Hamiltonian that describes it is time independent, hence the idea 
of using rotating coordinates in the present case cannot be applied. Moreover 
our problem is nondegenerate, however even in our case it is advantageous 
to perform a change of variables and apply a variation of the action angle 
variables used in Refs. 4, 12. Here the nondegeneracy of the problem plays 
a role similar to the rotating coordinate system in the restricted three body 
problem. 

For the 5*3 symmetric orbits we can instead consider the equations in the 
rotating frame, and prove a theorem similar to the one proved in Ref. 8 
for the anisotropic Kepler problem (in Ref. 8 the author remarks that the 
analysis of the Kepler problem can be redone in the Manev case, but he 
doesn't provide a proof) . 



III. The Si Symmetric Orbits with i = 1, 2 

We recall that the action variables introduced in Refs. 4, 12 are given by 



where h is the energy constant and K is the angular momentum. These 
variables are defined for /i < and K'^ > 26, 7 > 0, to avoid coUision orbits 




(5) 
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as well as circular orbits. The related frequencies are 

< 



K 



and 6 and are the angle variables associated to K and / respectively. 
The unperturbed Hamiltonian in the new variables can be written as 



Ho 



1 



Now we can consider new variables that are linear combination of the previous 
ones. They are defined by the following canonical transformation 



(6) 



L^K + I 
G = -I 

I = 9 
- 9 = 0-(l> 

Where / is the mean anomaly (where l{t) = ujL{t — ^o) ^-nd to is the time of 
pericenter passage), g is the longitude of periccnter as they are defined for 
the Manev problem in Ref. 13. Moreover also the action variables can be 
written in terms of the orbital elements of the Manev problem. If we set 



1 



2\h\ 

as in Refs. 4, 13 then 



and e = ^Jl - 2(K^ - 2b)\h\ 



,1/2 



l-(l-e^) 



2N1/2' 



and 



where a is the pseudo-semimajor axis, e is the pseudo-eccentricity, and the 
sign + (resp. -) holds for X > (resp. < 0). The conditions to avoid 
collision orbits and circular orbits, on which g becomes meaningless, can be 
written in terms of the orbital elements as a > and < e < 1. The new 
unperturbed Hamiltonian is 



Ho 



1 



2{-G+^{G + L)^-2by 



(7) 
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so the perturbed equations of motion become 



L ■■ 
G 

i = 

9 -- 



d{Ho+eW) 
dl 

_ d{Ho+eW) 

da 

d{HQ+eW) ^ 
dL 

d{Ho+eW) _ 
dG 



— e- 



dW 
dl 

aw 

dg 

,, , ^dW 
/, J- 



(8) 



where W is expressed in the new variables and 

G+L 



UJL = 
^G = 



(-G+^(G+L)2-26)3^(G+L)2-2fe 



G+L-a/(G+L)2-26 



(-G+v^(G+L)2-2fe)3^(G+L)2-26 



With these preparations, i.e. the introduction of the action angle variables 
d^), we are well on our way to estabilishing the following result: 

Theorem 1 Let'y{t) be an Si-symmetric periodic orbit of the Manev problem 
with i = 1,2. Let the period be t and set e = fi — 1 with e <^ 1. Then there 
exists a r-periodic solution of the Anisotropic Manev problem 7e(^) such that 

7e(t)= 7^+0(6). 

PROOF: Let's consider an 5*1 symmetric orbit of period r = 2Tim/k (m, k 
relatively prime integers). We remark that, since the equations of motion are 
autonomous, we can reduce to study the symmetric orbits that have either 
the pericenter or the apocenter on the positive x axis at t = 0. 
If at t = 0, e = 0, the pericenter of this orbit is on the positive x axis, and it 
is crossing the x axis perpendicularly, we have 



^(0) = and /(O) = 0. 



(9) 



Since the periodic orbit is 5*1 symmetric, by Proposition 2, at the half period 
one has 

g{T/2) = rrni /(r/2) = kn (10) 
that follows from the solution of fBI) for e = 0: 



L = const. 

I = UJlt 



G = const. 

9 = uJat 
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Now if, for e 7^ we consider only 5*1 symmetric solutions of (^), it follows 
from the implicit function theorem that if the functional determinant 



I dg/dL dg/dG ) ^ 



at 



t = t/2 







(12) 



(13) 



then ([10[) would be satisfied for e > . To compute the determinant we can 
by analyticity substitute (p!0| ) into (|1T]) to find out at the time t = t/2 that 

^"'^/^'^ ^0 (14) 



D 



G + J{G + L)2 - 26)7((G + L)2 - 26)3/2 



Thus the existence of Si symmetric periodic orbits of period r obtained from 
the T periodic Si symmetric solutions of the unperturbed problem, that at 
t = have the pericenter on the positive x axis, is readily established. 
On the other hand, if at t = 0, e = 0, the apocenter is on the positive x axis, 
and it is crossing the x axis perpendicularly, we have 



^(0) = tt/A and /(O) = -vr/A 



(15) 



where A = {ul — ujg)/ujl- By Proposition 1, at the half period we have 

^(r/2) = (m + 1/A)7r /(r/2) = (-1/A + A;)7r. (16) 

Instead of computing the functional determinant directly, in this case, it is 
easier to consider the new variables given by the relations. 



L 



G 



G 

1 = 1 + vr/Ao 
[ 9 = 9-tt/Xo 



(17) 



that define a family of canonical transformations parametrized by Xo{Lq, Gq). 
For each orbit choose a different transformation from the family ([T7|), where 
Ao = A is a fixed quantity defined by the value of the action variables along 
the periodic orbit under consideration. 

The equations (|16D, expressed in the new variables, are of the same form as 
in (p^. Thus the functional determinant, in the new variables, is exactly D, 
and the existence of the remaining 5'i-symmetric r-periodic orbits follows. 
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Now the proof for the S'2-symmetric orbits can be done along the same lines. 
Consider an S2 symmetric periodic orbit of period r = 2iTm/k. If at t = 0, 
e = the pericenter of the orbit is on the positive y axis and it is crossing 
the y axis perpendicularly, we have 

^(0) = 7r/2 and /(O) = (18) 

Since the periodic orbit is 5*2 symmetric one has, at the half period 

g{T/2) = mn + n/2 /(r/2) = /tvr. (19) 

Now we consider only S2 symmetric solutions of (§) for e 7^ again it follows 
from the implicit function theorem that if the determinant D computed at 
t = t/2 for e = is non zero then (0) would be satisfied for e > 0. It 
is trivial to see from ([I^) that D ^ 0, and hence we found 5*2 symmetric 
periodic orbits for the perturbed problem. 

For the S'2-symmetric orbits having the apocenter on the positive x axis at 



t = the canonical transformation (|T^ can be used. Again we find the 



same expression for the functional determinant and hence, by the implicit 
function theorem, the existence of the remaining S'2-symmetric periodic orbits 
is proved. 

It is interesting to remark that Theorem 1 and its proof can be easly extended 
to consider any S'j-symmetric perturbation with i = 1,2 and a very general 
class of nondegenerate integrable Hamiltonians, however such a generaliza- 
tion is trivial and not strictly related to the problem under consideration and 
hence it will not be discussed any further. 

We can also observe that for 6 = 0, i.e. for the Kepler problem, the determi- 
nant in (|14D is zero. Thus in the case of the Anisotropic Kepler Problem, the 
continuation theorem proved above cannot be applied, and the existence of 
symmetric periodic orbits of "second kind" (for weak anisotropies) remains 
unclear. On the other hand the continuation theorem that we prove in the 
next section (for the circular orbits) can be applied to the Anisotropic Kepler 
Problem^ and hence at least the existence of symmetric periodic orbits of the 
first kind is a well established fact. 



IV. The S3 Symmetric Orbits 

Again we can consider the Anisotropic Manev Problem taking the parameter 
/i close to 1. Let the flow $(t, (r, i"), /i) of the equation of motion (1). In this 
section we prove the following theorem: 
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Theorem 2 Let r^{t) be a S'^- symmetric periodic orbit of the Manev prob- 
lem, i.e. a circular one. Set e = — and let r be the period o/r°(t). Then 
there exists a r-periodic solution $(t, (r(e), r(e)), e) of the Anisotropic Manev 
problem such that <l>(t, (r(0), r(0)), 0) = (r°(t), r°(t)). 

A. The equation of motion 

Now using the same notation as in Ref. 8 let r°(t) be a circular solution of 
the Manev problem which correspond to /i = 1 in the xy-plane, uj its angular 
speed and a its radius. For e = — 1 7^ we set, 



Expanding VH in powers of /i — 1 sufficiently small, after substituting the 
expression for r given above, considering the notation r°(t) = x^it) + iyoit) 
and s = u + iv we have that r(t, e) is a solution of equation of motion defined 
by if, and only if, s(t, e) is a solution of the equations 



r(t,e) = r°(t) + e s(t,e; 



(20) 




^ + ^)t; + r/(t) + 0(e) 
¥ + i)^ + eW + 0(6) 



(21) 



aS ^ a<i 





COS tut + i sin tut 



62 = «ei 



and using the same notation as in Ref. 8 where 



s = xiei + X2e2, 



s = 



equation (|2l|) can be written in an equivalent form as: 



z = Ao{t) + Az + 0{e), 



(22) 



where z = (xi, ^2, 2/2)'^, and 



11 



^0 






a{t) 

V m J 



/ 





UJ 


1 


\ 




—cu 








1 












U! 


V 







—cu 


/ 



where 



or equivalently, 



a{t) cos cut — P{t) sin cut = rjit) 
a{t) sincut + (3{t) cosurt = ^{t) 



a{t) 

m 



sin2a;t(^ + § 



V2a2 ^ as) 

— sin Oi)t COS + |l) 



(23) 



The eigenvalues of A are 0, with multiphcity two, ^-nd — ^rr- One of the 
two eigenvalues vanishes because the system is autonomous, and the second 
due to the presence of the first integral H. 

Now consider the real Jordan form J of A. The matrix J is defined by the 
relation J = T~^AT where T is 







V 









a 

2 afi 





a 







-2 











and the columns of T are the generalized eigenvectors of A. 
The vector Jq = T~^Aq and the matrix J are: 



/ 3i{t) \ 

m 

hit) 
V hit) J 



J- 



/ \ 

10 

4 

V y 



3^2 _ t£(t£2a4+26)v_i 



where the fact that ji{t) = (2iu a — 
about Jo that we need to retain. Furthermore we remark that uj'^a'^ — a — 2h 



) ^P{t) is the only information 
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gives the relation between a and uj and solving this equations gives only 

one positive solution (for 6 > 0). 

Letting z = T(, the equation of motion becomes 

C = Mt) + JC + Oie), (24) 

and its flow is given by 

^(t,C,e) = 7W + e-'* + 0(e) (25) 
where by the variation of constants 

-f(t) = e-^' [\-^'Ms)ds (26) 
Jo 

Therefore we have 



and from ( p6l) 





/ 1 













\ 
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1 















COS 4t 




sin 
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sin 




COS^t 


J 


we obtain 






( Hit) 


\ 










73(t) 














I 74 (t) 


J 







(27) 



where we retain only the information that 

71 t = 2cuV - ^ T^^) / f3is)ds. 28 

Jo 

B. The periodicity equation 

Since the Hamiltonian H of the anisotropic Manev problem is S'3-symmetric, 
as we have shown, we can write the periodicity equation as in^, 

$(I (r,r),e) =-(r,f). (29) 
13 



Then it easy to check that (r, f), e) is a periodic solution of the equation 
of motion with period r. To find periodic solutions we have to verify that (pQ]) 
is satisfied for a family of initial conditions. Equation (^) in ( coordinates 
is 

(30) 



T 



where ip(t, (, e) is the flow of (pll). Let us denote by V{(, e) the left hand side 
of the periodicity equation (^91) , that is, let 



V{C, e) = ^ (r/2, C, e) - C = i{t/2) + {e'^ - l) 
Using (p5|) we notice that the requirement 

7(r/2) + (e-^^ - l) C = 0, 



0. 



v{c,o) = 

imposes the restrictions 

7i(r/2) = Cr 



--72 (r/2) C2 = arbitrary 



and 



C3 = 2(i-cos.>) (-73(r/2)(cosa* - 1) + 74(r/2) sin a*) 
C = 2(i-;oL*) (73(^/2) sina* + 74(r/2)(cosa* - 1)) 



(31) 



(32) 



(33) 



(34) 



where a* = 7r(l + 26/a)"i/2_ g^sy to see from (H) and (||) that 71 (r/2) = 
0, therefore, we take 

c = (cr,C2,C3,C4T, (35) 

with Q arbitrary, for the moment. Now using the flow , we determine 
that the Jacobian matrix of V with respect to the variables ( evaluated at 
the point {(*, 0) is given by 



/ \ 

r/2 

COS a* — 1 sin a* 

\ - sin a* cos a* - 1 J 



(36) 



Consider the system of three equations formed by those in (|3TD corresponding 
to the indices i=2,3,4 and fix the variable C2 = C2- Jacobian matrix has de- 
terminant r(l — cosa*), that is always positive since < 7r(l + 26/a)~^/^ < tt. 
Therefore the implicit function theorem guarantees the existence of analytic 
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functions Q = Cj(^); = 1,3,4 in a neighborhood of e = 0, satisfying the 
equations 

P,(C,e) = 0, (^ = 2,3,4) (37) 

where 

C(e) = (Ci(e),C2,C3(e),C4(e)) (38) 

and such that 

0(0) = C (^ = 1,2,3,4). (39) 

It remains to show, in order to have periodicity, that also the remaining 
equation 

Pi(C(e),z/(e),e) = 0, (40) 

is satisfied in a, possibly smaller neighborhood of e = 0. That will be done 
employing a first integral of the system under discussion, i.e. the Hamilto- 
nian. 



C. Integral of motion 

Since the Hamiltonian is a integral of motion of the problem under discussion 
we can apply the same analysis as in Refs. 8, 14. In particular using the 
same notations as in Ref. 8 we can define 

where H^{z,t) is a time- dependent, r-periodic first integral for system (P^. 
The above integral satisfies the following relation 

H,iz,t + T/2) = H,{z,t) (41) 

for all t, since H{—r, — r) = H{r, r), r(t) = r°(t) + e s{t) and 

r°(t + r/2) = -r°(t) , s(z, t + r/2) = -s(z, t). 

Performing a change of coordinates we can define 'HeiCi t) = H^{T(, t), hence 
(|4l| ) can be written as 

n,{C,t + T/2) = n,{C,t). (42) 

Moreover since Tie is an integral of motion it verifies that 

7^,(0(C,e,t)) = iJ,(C,O). (43) 



15 



Thus applying equations it follows that 

7^,(^(r/2,C,e),0) = 7^,(C,0) 

and by means of the Mean Value Theorem we obtain 

Vc7^e(C,0)-P(C,e) = 0, (44) 

where V^Tie is the gradient of Ti^ with respect to C, and C is a point on the 
segment joining ^ to ^/'(t/2, ^, e). 

Expanding \l/(e) = tp{T/2,(,e) in power of e sufficently small it is easy to 
show (see Ref. 8) that \E'(e) = C* + 0(e) and consequently 

C = <(e) + (l-s)M/(e)=C + 0(e) 

for some s G (0, 1). Moreover if we also expand the Hamiltonian H^{z,0) in 
powers of e we get 

H,{z, 0)=Ho + t{Hr + H2-z) + 0{e^) 
or, in ( coordinates 

H,(C, 0) = no + e{7ii +^2-0 + 0(6^), (45) 
where Ho = Hq = Uu'^a^ - i - A), = and = H2 = T^{a~^ + 



a a2 

2ba~^, 0, 0, aui) = {au'^Ci, 0, 0, 0). Hence we obtain 



^Vc^,(C,0)= 7^2 + 0(e). (46) 



With these preparations equation ([44f) reduces to the equation in the un- 
known Vi 

[au'^ + 0{e)]Vi = 0. (47) 

since, for small e, we already found in Sec. IV B that "Pj = for z = 2, 3, 4. 
It is easy to see that for e = the equation above has solution Vi = 0. 
Thus, by continuity, [007^ + 0(6)] is different from zero for e sufficiently small. 
Therefore for such values of e this equation has a unique solution that is the 
trivial one. Consequently the remaining equation 

^i(C(e),e) = 0, 

is also satisfied in a possibly smaller neighborhood of e = 0. Hence all the 
equations of the periodicity system (|3lD are satisfied when ^ = ^(e), as long 
as e is sufficiently small. This completes the proof of Theorem 2. 
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